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/aSTR^T 

In this thesis, the problem of finding the mean and 
standard deviation of the project duration in PERT networks 
is considered. The conventional method, which assumes the 
expected values of activity durations as actual times for 
completion, always gives lowest bound to the project duration. 
Monte Carlo Simulation technique, which can give fairly accu- 
rate results, takes very high time even for networks of 
reasonable size, A new approximate method is developed in 
this v/ork which appears to be a reasonable compromise between 
the two extremes. The method enxjmerates k-th longest paths 
in the network and finds the maximum of corresponding random 
variables. , Dependenc^r among the paths is also considered in 
an effective way. The suggested method is extensively tested on 
randomly generated PERT networks with various topological 
structures. The results are compared with conventional method.^ 
simulation method and Sculli's method. The results' of the com- 
putational analysis, v/hich are found to be quite encouraging, 
are also presented. 


CHAPTER I 


INTRODUCTION 

1*1 INTRODUCTION! 

One of the most important problems in the analysis of 
PERT networks is the determination of the total completion time 
of the project . When the completion times of activities of the 
project are random variables, the project duration is also a 
random variable with a distribution f motion that is a complex 
function of the corresponding distribution functions of each 
activity. Because of this complexity in evaluating the distri- 
bution fmction, the calculation of mean and standard deviation 
of the project duration has become a difficult problem and 
drawn the attention of many researchers. There is no exact 
method or formulae available in the literature for finding the 
mean and standard deviation of project duration and one has to 
content with approximate methods. 

The conventional method of solving the problem assumes 
the time of each activity as equal to its expected value and 
treats the network as a deterministic one. It is easy to realise 
that this gives only the crude lower bound on the project duration. 
Monte Carlo Simulation can be used to get fairly accurate results 
but it takes a very high time even for networks of reasonable 
size. Hence the research work has been, directed to develop 
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approximate methods which improve the results of conventional 
method in a reasonable amount of time. Some of the existing 
approximate methods are due to F-ulkerson [ 2 ] , Elmaghraby [4] , 
Robillard and Trahan [6], Sculli [s] etc. Though these methods 
are able to improve the accuracy to some extent, there is still 
a great scope for work to further improve the accuracy in the 
mean and standard deviation of project duration. 

This thesis work is directed to developing a new approxi- 
mate method for the problem which gives tolerably accurate results 
even at an acceptable increase in the computational time, 

1.2 PRELIMINARIES J 

The following terminology, notations and statistical 
rules are used quite often in the thesis work. For the sake of 
completeness, they are given below, 

1,2,1 Networ k Notati ons i 

A directed network G(N, A) consists of a finite set N 
of elements, called go das, and a set A of ordered pairs of nodes, 
called arcs_, with some weights attached to the nodes and arcs. 

Let, n = In] = no, of nodes in G 
m = jAl = no, of arcs in G 

If there is an arc e from node i to node j, then i is 
called Tai l node, of the arc e and j is called Head node . Also, 
the arc e is said to be i ncid ent from node i and jfi cident to 
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node Let and 0^ denote, respectively, the sets of arcs 
incident to and incident from node i. The node S e N which 
has no arcs incident to it, is called the Source and the node 
t £ N which has no arcs incident from it, is called the Sink , 

A Path in G = (N, A) is a sequence i^,i 2 ,. . .,ij, of 

distinct nodes of N such that either (ij^, £ A or 

(ik^l» ijj) ^ A each k *= 1,2,..., r-1, A directed Path is 
defined similarly, except that (i^, a A for each 

k =» 1,2,..., r-1. A cycle is a path together with an arc 
(i^, i^) or (ij, i^). A di rected cycle is a directed path 
together with the arc (ij,, i^). . 

A network is said to he acyclic if it does not contain 
any cycle. Topological Sorting is a systematic way of finding 
the existence of a cycle in the given network. The procedure 
for topological sorting is as follows l 

Step 0: Set i = 0 

^KArmnmmm, uihr-m n,m\m m 

Step 18 Set i = i+1 

Identify a node with no incoming arcs and number it as i. 
Delete the node i and its outgoing arcs. 

If i = 'f then STOP, else go to Step 1. 

It is obvious from the. above procedTore that all nodes 
can be numbered if and only if there is no directed cycle in 
the network. 
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A directed complete network is one in which for every 
node i £ N, (i, k) £ A for all k = i+1,.,,, n. The Length of 
a path in G is the sum of the lengths of all the arcs on that 
path. 

1.2.2 PERT Terminologiri 

Any project consists of a number of job operations or 
tasks to be performed which are called activities . An Event is 
a specific instant of time which marks the beginning and the 
ending of an activity. 

Network representation of a project is very useful in 
Project I'fenagement because it gives a clear picture of the 
project in a ^apshot. In this representation, the activities 
of the project correspond to the arcs, the events to the nodes 
and the activity durations are the weights attached to the arcs. 
The position of an activity is dictated by the following rules 
all the activities associated with the incoming arcs to a node 
must be terminated before those associated v/ith the outgoing 
arcs at the node can start. In consequence, a project network 
can be only directed acyclic. 

If the durations of all the activities can be predicted 
exactly before starting the project, then the network is said to 
^ deterministic network . On the other hand, if the activity 
durations cannot be predicted exactly due to some chance varia- 
tions, then the netvrork becomes a probab i listic networ k or 
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PERT (Program Evaluation and Review Technique), This is often 
the case of research and development activities, projects of 
non-repetitive nature such as defence projects. 

Pr o j ect duration is the total time required to complete 
the project. The sequence of activities which decide the pro- 
ject duration form the critical path. The path for which the 
length is nearly equal to the critical path length is called a 
near_ critical pat h. 

1.2.3 Probahility Theory in PERT? 

The probability theory and mathematical statistics have 
found a great use in the analysis of PERT networks. One of 
their most important uses in PERT network analysis is to get 
a measure of the uncertainty that the project will be completed 
before a specified due date. The statements of the possible 
range of times for the completion of the project and the probabi- 
lities associated with each are quite useful for the management. 
By adding to this information an appraisal of the consequences 
of not meeting a scheduled date, the management can better plan 
a project. 

Normal distribution is often referred in the thesis 
work. It is a symmetrical bell shaped curve with two parameters, 
mean (m) and standard deviation («"), The probability density 
function for a normally distributed random variable X is 
given by. 
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f(x) = 


e 


1 

2 


(^) 






The cumulative distribution function, F(x), involves the inte- 
gration of the above function, between the limits “<» and x, 
which is difficult to evaluate. 


The normal distribution with the parameters m = 0 and 
<3- = 1 is called S tand ar d Normal distribution . Its density- 
function and distribution function, denoted by ^ (x) and 0(x) 
respectively, are given by, 


v(x) 



e 


1 2 
2 


$ 


and 0(x) 


1 2 
X “ x^ 

r e dx . 

^ -oo 


The importance of 0(x) is that it gives the cumulative area 
upto the point x in the Standard Normal Curve, i.e, it gives 
the probability for an observed value of X to be less than or 
equal to x. The values of 0(x) are tabiolated for different 
values of x and are made available in any statistical book. 


As a general case, if one wants to compute Prob [x <_ x] 
where X is Normally distributed with parameters m and o-', then 
X can be converted to the Standard Normal Scale by x = (x“in)/o- 
and the tables for that value of x can be seen. It is interest 
ing to note that 99.7 percent of theobserved values for X lie 
within the range (m ~ 3<r') and (m + 3^^). 
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With the background about the Normal distribution, we 
now present some formulae and theorems available in the litera- 
ture. 

( ^ ) S l_^rkls , N ormu^^ [ l3] i 

These are used for finding the mean and standard devia- 
tion of the maximiJin of two independent Normal random variables.. 
Suppose and are two such variables with parameters 
and (m 2 , «~ 2 ^ respectively. If Y = max (X^, X 2 ), the parameters 
of Y, m and <r , are given by 

m = m^ 0(a) m 2 0(-oc) + ay (a) 

V2 = (m^^ 0 (a) + (m2^ ^^ 2 ^ 0 C*-cx) 

+ (m^ m 2 ) a "V'(a) 

2 

and, cr- = ■V 2 - m. 

where, a = '*'^2 “ ” ^^1 ^ 


( ii ) A pproximat e. F ormul a fo r 0 (a ) ■ 


The following approximate formula given in Greer and 
Cacava [ 9 ] can be used to find 0(x) for a given value of x. 


0(x) ~ 1 


(l+c-,x)^ (l+c-,x)^ (l+c^x)^ 


(l+c^x)^ (l+c^^x) 

where, = 0,2316419, = 1.330274, c^ = 1.821256 


c^ = 1.781478 , c^ = . 0,3565638, Cg 


0.3193815. 
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It is important to note that if x is negative, the 
above formula should not be used directly. Instead, it can 
be used for the absolute value of x and then use 0(-x) = 1-0 (x). 
This formula gives 0(x) upto 5 digits of accuracy. 

(iii) Jensen's Inequality [l5]i 

According to this inequality, the expected value of the 
maximum of a set of independent random variables is greater 
than or equal to the maximum of the eixpected values. Mathe- 
matically, if X^, X^, ..., X^ are independent random variables, 

E[Max (X^, X^, ..., X^)] > Max [E(X^), ECX^), ..., E(X^)] 

Similarly, the variance of the maximum is, in general 
less than the minimum of the variances. This is particularly 
true when the individual distributions are symmetric such as 
Normal. 

( iv ) Centra l_ Li mit T h eorems 

This theorem is one of the most in^jortant theorems in 
all of mathematical statistics. According to this theorem, the 
distribution of the sum of n independent random variables is 
approximately Normal regardless of the distributions of indivi- 
dual random variables, if n is sufficiently large, 

(v) Sums of Random Var iables s 

The mean of the sum of n independent Normal random vari- 
ables is equal to the sirni of the means of individual random 
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variables. Similarly, the variance of the sum is the sum of 
the variances. 

Mathematically E [Z X^] = EE (Xj^) 

Var [£ X^] = E Var (X^) 

1.3 MOTIVATION OF THE- PROBLEMS 

The problem of finding the project duration is very easy 
in deterministic netv/orks as all the activity durations are 
fixed and well known, Infact, the project duration becomes the 
length of the longest path in the network. 

The same problem, hov/ever, becomes exceedingly^' difficult 
in PERT networks. Since the activity durations are random, vari- 
ables, each path in the network from source to sink also becomes 
a random variable which is equal to the sum of all the random 
variables corresponding to the activities on that path. Hence 
project duration also becomes a random variable that is equal to 
the maximum of all the random variables corresponding to the 
paths. The mean and the standard deviation of project duration 
are difficult to evaluate because of two main reasons (i) the 
number of paths in the network may be intractable for large net- 
works and (ii) we need to find the maximum of a set of dependent 
random variables, corresponding to the paths in the network, for 
which no method is available. 
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The conventional method assumes the expected value of 
the activity duration as the length of the arc in PERT network 
and treats the network as a deterministic one. The mean and 
variance of the longest path in the network are treated as the 
parameters of project duration. If there are more than one 
longest path in the network, then it takes the path with highest 
variance as the longest path. Thus, the conventional method 
assumes the maximum of means of various paths as the mean of 
Project duration. But hy Jensen's inequality, this is only a 
lower bound. Thus conventional PERT method is not acceptable 
when the accuracy in the mean and Standard deviation of Project 
duration is required. 

The importance in the accuracy of mean and Standard 
deviation of Project duration comes from the fact that these 
parameters are used in computing the probabilities of meeting 
some specified due dates and if these parameters are inaccurate, 
the probability statements will be very much misleading. To 
get a feeling of the errors in the probability statements, let 
us consider the following simple example 

Consider a project for which the mean and standard 
deviation of Project duration are found to be 22.00 and 2.00 
respectively by conventional method. Then the probability of 
completing the project in 21 uinits of times is = 0.5000, 

Suppose the actual mean and standard deviation are 23.00 
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and 1,50 respectively. Then the probahility of this event is 
0 _ 0.0918. In other words, v/hen there are only 9.18 

percent of chances for completing the project in 21 units, the 
conventional method indicates it as 50 percent, ¥e thus see 
that the statements can be highly misleading even when there are 
only small errors, in the parameters of project duration. 

The errors in the parameters of project duration 
become higher if there exist near critical paths in the network. 
For example, consider a network for which the distributions of 
various paths are as shown in Fig. 1,1, The critical path 
has parameters (m^, Consider a near critical path P2 

with parameters (m2,"~2^* Now, if bad luck is generally experi™ 
enced for the activities along P2 and good luck for the activities 
along P^, it is quite possible for P2 to exceed P^ in length — 
that is to become critical itself. Infact, the same may happen 
with other near critical paths too. Therefore, it is overly 
optimistic to assume (m^^, as the parameters of Project dura- 

tion as assumed in the conventional PERT. 

The main idea in the present work is to reduce the errors 
in the parameters of Project duration by considering the effect 
of near critical paths. It is assumed, without loss of genera- 
lity, that all the activities in the network are independent. 

It is also assumed that the durations of individual activities 
are Normally distributed. This assumption has often been made 
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in the literature and this can be justified by the fact that 
most of the networks can be reduced to a guide network, where 
a completely independent path becomes one activity. Then the 
Central limit theorem justifies the Normality assumption for 
the activity durations in the guide network. 

1'.4 FORMULATION OF THE PROBLEMS 

¥e nov;- give the Statistical formulation of the problem 
which forms the basis for our further discussions. Let us 
denote the various terms as follows t; 


hd 


i random variable corresponding to the 
duration of activity (i, j). 


(m. V. .)i Parameters of t... 

-^0 -*-U 

P, 


k 

Pk 


s Set of all the arcs on k-th path, 
s Random variable corresponding to P^» 

(Mk, s Parameters of p^. 

1 s Random variable representing the project 

duration. 

Now, Pj^ is given by, 

p,. = I t.. 

(l,d)eP„ ^ 


therefore, = E(p^) = E ^ 


kk 


I E(t,J 


(i,j)ePj^ (l,d)ePv 


Z m . . 


k 



14 


Similarly, 


V 


k 


S V. . 

(l.d)eP^ 


The Project duration 1 is given b 3 ^ 


1 = Max (P^, P^j Pp} . 

We are interested in approximating the parameters of 1. 
r.5 OUTLINE OF THE TI-ISSISs 


A brief, chapter by chapter, outline of the thesis is 
given belowa 


In Chapter II, the existing literature related to the 
problem is discussed. The following methods are described in 
some detail. 

(i) Malcolm's algorithm 

(ii) Fulkerson's algorithm 

(iii) Elmaghraby's algorithm 

(iv) Martin's algorithm 

(v) Robillard and Trahan's algorithm 

(vi) Sculli's algorithm, and 

(vii) Simulation approach. 

Chapter III is concerned v/ith the new approximate 
method developed for the problem. The intuitive ideas which 
underline the development of the method are first described. 
Then a stepwise description of the algorithm is presented. 
Lastly, the method is illustrated v/ith a numerical example. 
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In Chapter IV, the computational investigations of the 
suggested method, together with the conventional method, simu- 
lation approach and Sculli’s approach, are presented. The 
results are presented for a large n-umber of network problems 
generated randomly. 

A list of references and bibliography is given at the 
end, which is follov/ed by program listings of all the programs 
used in the computational study. 



CHAPTER II 


LITERATURE SURVEY 

2.1 lOTRODUCTIONi 

During past several years, the problem of finding the 
project duration in PERT networks has received widespread 
interest. Most of the research was directed towards obtaining 
an approximate estimate because of the inherent complexity 
invo3.ved in evaluating its exact density function. Some 
researchers concentrated on finding only the mean of the 
project duration while others on finding both the mean and 
standard deviation. 

In this chapter, we briefly review various approaches 
available in the literature. The gradual improvements in 
each method are discussed and their limitations are pointed 
out . 

2.2 MALCOLM'S APPROACH! 

The first approximation for the expected completion 
time of the PERT network was proposed by Malcolm [l] , who 
-suggested to use the expected values of the individual activi- 
ties in evaluating the length of the longest path in the network 
The parameters of the longest path are considered as the mean 
and standard deviation of the project duration. Hence the 
expected value of the project duration, g^, is evaluated by 
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defining a function recursively as follows s 

gl = 0 

Si = [Sfc + 

where is the set of incoming arcs at node i, and 

t^i is the random variable corresponding to the arc (k, i) . 

Thus the model is reduced to a deterministic ^orm with- 
out incorporating the Stoohastic element, variance of activity 
duration. That is v/hy, in many situations , this estimate is 
observed to be far from the real value [l6] . 


This approach is referred to as the conv entional. , meth od 
because of similar calculations involved as in deterministic 
networks . 


2-. 3 FULKERSON'S APPROACH? 

In an effort to improve the above estimate, Fulkerson [2] 
proposed, for discretely distributed random variables t. ., the 

J. J 

function f obtained recursively from 


% - 0 


f.. = 


E_ P(I^. ) max {fj^ + t^^} 


kelj^ 


^ki- 


where is the joint probability distribution function of 

all t, - , k e I . . 
ki’ 1 


Thus he treated t^^ as random variables in evaluating 
f ^ which allows a greater choice of path taken to have a length 
approximating the critical path length. Obviously, f is 
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uniformly better than and in particular 

< 1 

XI ■* n " n 

. where 1^ is the actual expected value of project dtoration, 

Fulkerson 's approach bears straight forward generaliza- 
tion to the case of continuously distributed random variables. 
It is a well Icnown fact that the cumulative distribution func- 
tion ox the maximum of a finite set of independent random 
variables is equal to the product of the individual cumulative 
functions of the variables in the set. This fact was utilized 

dy Clingen [3] to obtain a computationally tractable form for 
f as follows • 


h i-i 

f^ = b - J 


a k=l 


where 


a = max (f, + a, . ) 

kel. ^ 

1 

b = max (f^ + bj^. ) 


kel. 


and 


z-fk 




'ki 


where, p (x) = density function of t, . 

ki 

^ki "" minimum possible value for and 

\i maximum possible value for t, . 

Kll 

The limitation oi the Fulkerson 's approach is obvious 
from the expression for If the integral of the 
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density function of the distribution of any activity is diffi- 
cult to evaluate as in the case of Normal distributions, then 
this approach cannot be used. Also, this approach assumes all 
fj^, k e as independent. But it is easy to see that they 
may be dependent because of the common arcs in various sub-paths 
reaching to nodes k e 1^^, 

2'. 4 ELMAGHRABY'S METHODS 

Elmaghraby [4] extended the Fulkerson's approach in 
order to obtain still better approximation. His approach is 
based on the following simple observations if all the arcs 
in a directed acyclic network are reversed, the average project 
duration remains unchanged. However, intermediate values of 
f^, 1 < i :< N, donot necessarily remain the same. Consequently, 
f^ can only bo improved if we take the maximum of the two values 
obtained from 'as given' and 'reversed' subnetwork. 


Let be the estimate using Fulkerson's approach in 
the reversed subnetwork. Let s^ and Uj|^ be defined recursively as 

= 0 

P(I^) max {sj^ + t^^} 


Si = U3_ 


Ui = 


kel^ 


s^ = max (u^, 

s^ is the estimate of 1^ and obviously 





20 


2.5 MARTIN’S APPRO ACHi 

Martin [l7] presented a computational method for the 
evaluation of the density function of the project duration in 
a PERT network under the assumption thcit the arc duration 
density functions are polynomials. This approach consists of 
mainly two steps as follows s 

The first step is to transform any directed acyclic 
network into a series parallel netv/ork. Two subnetworks are 
said to be in series if the Sink of the first one is the Source 
of the second one. They are in parallel if both have same 
Source and Sink. This step is performed by depicting the PERT 
network as a tree where some nodes and arcs are duplicated as 
necessary. 

The second step is to reduce systematically a series 
parallel network to a single equivalent arc with a density 
function that characterizes the project duration, 

Martin developed algorithms for performing the above 
two steps. Unfortunately, the algorithms require a great deal 
of calculations even with the assumption of polynomial density 
functions for activit 3 ^ durations. 

2.6 ROBILLARD AND TRAHAN'S APPRO ACHs 

Robillard and Trahan [6] generalized the Fulkerson's 
estimate by ejctending the concept of I^ to include all the 
subpaths in G terminating at node i. Let be the set of 
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subpaths terminating at node i and t^ be a random variable 
associated with a subpath p representing the total duration of 
the set of activities in p. Then 

= 0 

m. = max {m, + t } 
peP^ P 

where k is the beginning node of p. 


It is obvious that f^ general, it 

is difficult to evaluate m. because we have to handle many t 

1 P 

simultaneously. So they proposed two modifications in the 
above approach which can deal only v/ith arcs. Both these 
modifications require to reduce the subnetwork till node i to 
a parallel network as done by Martin. 


The first modification is based on the observation 
that if we permute the arcs along a path in the reduced PERT 
network, the average duration of the partial project remains 
imchanged. However, if v/e calculate f^ for various permutations, 
it does not necessarily remain the same. If the maximum value 
of f^ over all the relevant permutations is chosen, it can 
only be better than Say this estimate as n^. Therefore, 


The second modification is to reverse all the arcs in 
the reduced network, as done by Elmaghraby, using n^ instead of 
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If is the hound obtained by this procedure. 


^i ®i ^ V 


Hi- 


The above approaches require a lot of time for convert- 
ing the subnetv/orks into parallel networks and for checking 
various permutations of arcs. 


Robillard 'and Trahan [?] have also suggested another 
method for evaluating the distribution function of project 
duration using the approximation proposed’ by Davies [l8] in 
calculating the inverse characteristic function. 


2.7 SCULLI'S METHOD I 


As it is already discussed, Fulkerson 's approach cannot 
be used when the distributions of activity durations are Normal. 
Sculli [s] developed a method to overcome this limitation of 
Fulkerson's approach. He asstimed that the durations of indivi- 
dual activities are Normally distributed and used Clark 's 
approximate formulae for finding f. - max {fi- + 'tv,-}* 

It is assumed here also, as in Fulkex'son's approach, 
that all f^, k £ I^ are independent. As a result, the errors 
in fj^, particularly in dense networks, are more because of 
large number of dependent paths terminating at nodes k e 1^^. 
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Many other methods based on similar ideas have been 
proposed. Lindsey [5] developed an estimate for the project 
duration based on a model approximating the actual network. 

This estimate may fall on either side of actual value 1^, Its 

calculation is more complicated and it did not improve the 

estimate 

Fulkerson 's /significantly except for nodes nearer to Source 
node. Kleindorfer [20j tried to find the distribution functions 
that bound the activity starting and completion. In practice, 
this approach concerns only PERT net\TOrks where the distribution 
functions associated with the activities are discrete. 

2.8 SIIiULATION APPROACHs- 

A very close approximations to the mean and standard 
deviation of project duration in PERT networks can be obtained 
by using Monte Carlo Simulation [l9]. Given the distributions 
of activity durations, the network can be simulated by generat- 
ing the durations of activities randomly and finding the project 
duration using conventional method. In each run, we get a value 
for the project duration and the mean and standard deviation of 
all the values obtained for different runs gives the parameters 
of the project duration. 

The advantage of this approach is that it considers all 
the possible variabilities in activity times if the network is 
simulated for sufficient number of runs. So this method can 
be taken as an accurate method for finding the parameters of 
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project diiration. The major drawback in this method is 
that it requires a lot of time for simulating even a reason- 
ably large network. 



CEAPTER III 


PATH ENUMERATION ALGORITM 

3.1 INTRODUCTION! 

In this chapter, a new approximate method for finding 
■ the mean and standard deviation of project duration in PERT 
networks is developed. The. method is called Path Enumeration 
Algorithm (PEA) as it is "based on the enumeration of k-th 
longest paths in the network. ¥e first describe the intuitive 
ideas which led to its development. A stepwise description of 
the algorithm is then presented. A numerical example is also 
given to illustrate the algorithm, 

3.2 DEVEL0P]>1ENT OF PATH ENUMERATION ALGORITms 

As already discussed in Section 1,4, the actual project 

duration in PERT networks is given by 1 -- Max {p -} where p - is 

i 

the random variable corresponding to the i-th path in the network 
which is assumed to be normally distributed with parameters 
(Mi, V^) . The random variables p^^’ s are dependent because of 
the common are^s in different paths in most of the networks, ■ 
There is no method available to find the parameters of maximum 
of a set of normal random variables even for independent variables. 
Dependency among the paths makes the problem more difficult. Also 
the number of paths in a network may be extremely large even for 
networks of reasonable size. Hence it is quite unlikely to 
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enumerate all the paths in the network. It is, infact, not 
necessary to enumerate all the paths because only a set of paths 
generally dominate all other paths in the network. In the path 
enumeration algorithm, we enumerate this set of paths and 
heuristically try to find their maximum considering the depen- 
dency among themselves. 

It is easy to see that the paths v/ith higher value of 
expected path length tend to dominate those paths with lower 
values of expected path lengths. We therefore enumerate the 
paths in the non- increasing order of their expected lengths and 
consider paths until further enumeration changes the parameters 
of maximum of such paths by negligible amount. 

The enumeration of paths in the required order can be 
efficiently done by Dreyfus method [l4] , A procedure for find- 
ing the parameters of maximum of the random variables corres- 
ponding to the paths, v/hich are Normally distributed is yet to 
be discussed. One obvious method is to use Clark's formulae 
recursively. Simulation technique can also be used for this 
problem. In each simulation run, we can find the maximum value 
of the variates generated randomly v/ith the parameters of the 
random variables and then the mean and standard deviation of all 
such values for various runs gives the parameters of maximum of 
the random variables. Both the above methods were programmed 
and the values of the parameters were found to bo comparable. 
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Hence Clark's method, which takes lesser time than simulation 
is selected to find the parameters of maximum of random vari“ 
ables. 

Thus the procedure for finding the parameters of project 
duration is to enirnierate k-th longest path in the network and 
to find 'the parameters of maximum of this path with the previous 
paths using Clark's formulae and repeat these two steps until 
the enumeration of a path does not change the parameters signi- 
f icantly. 

The results from the above procedure were found to be 
substantially different from those with simulation approach. 

These deviations are anticipated since Clark's formulae treats the 
random variables as independent in finding the parameters of 
their maximum. But these random variables, which correspond 
to different paths in the network, may be dependent as there 
may be some common arcs. Until v/e consider the dependency among 
the paths in an effective way, accurate results may not be 
obtained, ¥e now describe a method to take care of dependency 
among the paths. 

Any two paths, and P^, in a network may have some 
common arcs and some uncommon arcs. If the sets of uncommon 
arcs are separated, then the two paths may be as shown in. Fig, 3.1. 
To consider dependency between and P^, the parameters of the 
naximum of only the uncommon portions can be found and added 
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to the parameters of common portion. We describe this in a 
more formal manner. Let, 



S 

= Set 

of 

arcs in P^ but 

not in 

h 

1 — 1 

II 

- h 

n 



P 

= Set 

of 

arcs in P 2 but 

not in 

h 

= h 

H 

1 

n 

h 

and 

I 

= Set 

of 

common arcs in 

P^ and 

h 

= h 

n p^. 




The arcs in the set I can be easily found by labelling 
all the arcs on the first path P^. The parameters of the common 
portion can then be calculated as. 


and. 


D, 


I ra. . 


Z V. . 
(i,a)el 


If ^ and n denote the random variables corresponding to the 
uncommon portions, their parameters are given by, 


(Cl. Dp = (Ml - C 3 , Vi - Dj) 

and (C^, D 2 ^ = ^^"^2 “ ^3’ ^2 “ ^3^’ respectively. 

Now the procedure for finding the parameters of maximum 
of tv/o dependent random variables corresponding to the paths 

and is to find the parameters of maximum of the independent 
random variables C and H using Clark's formulae and add them to 
(C^, D^). 

This procedure is to be extended for finding the maximum 
of a set of dependent random variables in. order to make it usable 
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in PERT networks. The difficulty in extending it is that a 
nev/ path may have different arcs in common with different paths 
enumerated and so it is difficult to decide which arcs are to 
he included in the set I, ha analogy v/ith the procedure for 
two paths in the network suggests to include the arcs which are 
common to the new. path and to any of its previous paths. So all 
the arcs on the previous paths are to be included first in the 
set S and then the arcs in set I can be easily identified by 
labelling all the arcs in the set S, 


The effect of the arcs in the set I, which are causing 
for the dependency, is now to be considered. The enumeration of 
a new path P changes the pare^meters of the maximum of all the 
previous paths (M, ) to certain values say (M^, V^) with 

Mr > and Vr < V^, Thus the affect of the arcs in P is to 
increase to Mr and reduce to Vr. In order to consider 
this effect in the further enumeration, the mean of these arcs 
is to be reduced by a factor of M^/Mr and the variance to be 
increased by a factor of call the calculated 

means and variances using these factors as effective me ans, 
and eff ecti ve vari a nc es, v^^ respectively. Therefore, 


m . . = m . . X 

10 10 






7 X vyvg 


^ (i,o) e 


¥ (i,0) £ P. 


The above formula.e require the initialisation of m. . and v. . 
for each arc in the netv/ork as equal to its actual mean and 
actual variance. 
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The effective means and effective variances are now 
to be used in finding the parameters of the common and uncommon 
po rt ions . Thus , 



Z 

(i, d)£T 



and 



Z 



~ and 

In the above expression for the term signifies 
the variance of all the arcs in the set S and signifies the 

variance of the labelled arcs in the path P. is calculated 

as the variance of maximum of all the previous paths but is 
calculated, in a different way, as the sum of v\j for the 
labelled arcs. Thus D_ is not in the same scale of V^. To 
bring it into the same scale of V^, we use a scaling factor of 
where is the sum of v^. for all the arcs in S. Thus, 


D- 


where , V, 


D, X (VVVj 


1/ ’s" » 

V . . 


(i,j)eS 

the scaling factor is less than or equal to 1. 

We thus have a heuristic procedure for finding the para- 
meters of maximum of random variables, taking care of dependency 
among the paths. This procedure can be repeated, after the 
enumeration of each path, until the difference between the 
successive values or the parameters of project duration is less 
than a small number, e, for a number of times. 
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3*5 STATEMEM' OF THE ALGORITHI'T? 

Path Enumeration Algorithm can be described stepwise 
as follows? 





Enumerate first longest path P using Dijkstra's 


algorithm. 

Let, 




ill = 

L m. . 

(i,j)£P IS 

and 

=r 

2 V. . 

(i,d)eP 

Set, L. . 

10 

= 1 V (i 

,d) e P 



m. . 

10 

= m. . ¥ (i^ 

,o) £ -A 



V. . 

10 

= V. . (1, 

I i ) £ A 



V 

s 

= \ , s = 

= P and 

k 

= 2. 

Enumerate the k-th longest path P using Dreyfus method 
Let, 

X = 

L m. . 

(i,j)£P 




If (X = oo) ^ 

then go to Step 9, else 

go 

to Step 3, 

Bef ine , I 

= s n P 




Set, = 

\ and 




H 

y mV- 

(i, 9 )eP mj 

and V 

= 

,2 V. . 

(i,3)£P 

■Define the parameters of 

the common 

portion I as 

follov/ss 





S = 

r 2 ra' ■ 

Ci,a)el 

and 



D-, = 

2 V. . : 

(i,o)£l 

X V^/V 

1' s 
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Define the parameters of uncommon portions as follov;s« 


= - C. 




- M - 


Dp = \ 


D, 


V 


D. 


If =: 0 then go to Step 8, else go to Step 4. 

Find the parameters (C^, D^) of the maximum of 
uncommon portions as follows® 

Set a - ^ and a = (C^-C^Va 

0(cc) + 0(“a) + a v(a) 

De = [(C^^ -i- D^) 0(a) + v D^) 0(-a) 

-h (C^ -h C„) a y{a)l - c/. 


Set I'h = C + C, 
1 e 3 


a ' ^2 

and V. 


1 


D + D_ 
e 3 


If (M^ “ > a) then go to Step 6, else go to Step 5 

Set count = count -i- 1 

If (Count >5) then go to Step 9, else go to Step 7. 


6i Set count 


0, 


Jji Update m. . and v. . as follows' 

J-J ij 


“ij = "■« V’t 


V (i,i) e p 


V... = . X V^/V^ ¥ (i,i) e P 


ij 


Update the labels and V as follov/ss 
L-. =1 ^ (i,j) £ P - P n s. 

X J 

V = V + Z V. . 

(i,d)£P“pns 


Set s 


SUP 
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step 8 s Set k = k + 1 

If (k >, Limit) then go to Step 9, else go to Step 2. 

Step^^^s and are the mean and variance of the project 
duration. STOP. 

3.4 NUMERICAL SXAi-IPLEs , 

We nov^ solve a numerical example to illustrate the 
various steps of the algorithm!. Since the main obstacle in 
computing the parameters of project duration stems from the 
possibility of non-critical paths becoming critical, the network 
shown in Fig, 3.2 is anal 3 /sed \'/here all the nine possible paths 
are expected to be critical. In this network, it is assumed 
that the individual activities have a Normal distribution X'/ith 
variance equal to 20 percent of their mean. The mean values are 
mentioned over the arcs in the network. 

The sets P and S at each iteration and the parameters 
of project duration considering the affect of all the previous 
paths upto that iteration are shown in Table 3.1. The iteration 
number indicates the number of longest paths that have been 
enumerated till that iteration. For this problem, the value 
of e is fixed as 0.2 percent of the mean of first longest path 
and the nximber for which the difference between successive values 
is less than s is fixed as 5. It is terminated after evaluating 
all the paths in the netv/ork because the condition with e is not 
satisfied. The parameters obtained at the terminating point;g 
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are Mean = 24.54 and Standard deviation = 1.43. The 

parameters obtained by Simulation approach for a total of 10000 
runs are 24,6? and 1.47 respectively. 



CHA TER IV 


COMPUTATIONAL INVESTIGATIONS 


4.1 INTRODUCTION! 

In this chapter, we present the computational performance 
of Path Enumeration Algorithm developed in Chapter III on randomly 
generated network problems. The resiUts obtained by simulation 
approach are considered as correct results as no exact method 
is available. 

The computational esqjeriments require a large number of 
randomly generated networks. The details of the network require- 
ments and the type of network problems considered for computa- 
tional study are discussed in Section 2. Section 3 covers the 
explanation of various algorithms used as subroutines while 
programming Path enumeration algorithm and other methods. The 
computational considerations such as computational times, accu- 
racy, storage requirements etc. are also given, for analysing 
the algorithms, in this section. In Section 4, the computational 
experiments with various parameters are explained and the 
resixLts of Path enumeration' algorithm are compared against 
conventional, simulation and Sculli's methods. ■ The advantages 
and limitations of Path enumeration algorithm are also given 
in this section. 
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4.2 TOPOLOGICAL STRUCTURES OF RAIDOM IJETWORKSs 

It is mentioned in Section 1,2,2 that the PERT networks 
can he only directed acyclic networks. Keeping this point in 
view, the following three random network generators are used in 
the computational study. These generate topologically sorted 
networks . 

^ ^ ^ Netw ork^ 

NG 1 takes the network dimensions, width (w) and 
length (li), as input and generates a sparse networks. An 
example network generated from NG l with a width of 3 and a 
length of 5 is shown in Fig. ,:4,1. There are 17 nodes and 
36 arcs in this network. Infact, one can derive the general 
formulae for finding the number of nodes, n, and the ntmiber of 
arcs, m, in the network as follows s 

n = ¥.L + 2 

m = L (3W - 2) + 1 

(ii) Networ k Gene rator 2 (N G Z)% 

NG 2 takes n as the input and generates a complete 
acyclic network. Therefore the networks generated by NG 2 are 
dense with m - C^) • An example network generated by NG 2 with 
n =* 6 is sho>/n in Fig. 4.2. 

(iii) Ne twork Generato r 3 (NG 3) a 

In the previous network generators, there was no control 
over the number of arcs. In NG 1, for a given width and length. 





both n and m are fixed. Similarly in NG 2, for a given n, m 
is fixed. In NG 3, we have control over both the parameters 
n and m. The number of arcs can also be specified in addition 
to the width and length of the network. Then NG 3 first gene- 
rates a skeleton network with the given width and length, and 
then it adds arcs between tv/o randomly generated nodes until 
the total number of arcs becomes the specified m. An example 
for the skeleton network is shown in Fig. 4,3 for a width of 
2 and length of 5, The final network, generated by NG 3, after 
augmenting the additional arcs for a given m = 1? is shown in 
Fig. 4,4, 

The computational study has been made with a large num- 
ber of problems generated by these three network generates . ¥e 
feel that these three generators are sufficient for our study 
and the results thus obtained will hold good for real-life- 
network problems . 

In all these network generators, the networks are first 
generated in Arc list representation and then converted into 
Adjacency array representation. In Arc list representation, 
the head and tail nodes of each arc, together with the arc 
lumber, are maintained. Therefore if there are m arcs in the 
letwork, the total storage requirement is 3 m. If all the 
ires incident from or incident to a node are to be retrieved, 
:he whole list of arcs is to be scanned. In Adjacency 
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representation, two arraj^s P0I3S1T and LIST are maintained. 

LIST contains all the arc numbers incidenting from all the 
nodes with the nodes arranged in sequential order. The position 
at which the arcs incident from a node starts, is stored in 
POINT array. The difference betvmen POIl'W (i+l) and POINT (i) 
gives the number of arcs emanating from the node i. Thus, by 
knowing the starting position in the LIST and the number of 
arcs, the arcs incident from that node are retrieved. Simi- 
larly the arcs incident to a node can be retrieved easily by 
maintaining two arrays RPOINT and RLIST. 


4.3 SOME COMPUTATIONAL CONS DERATIONS I 
^ ^ ) ^-gfl^^Q^tion of N o rmal, ra ndom, yariat es ° 


There are three well knov/n approaches for generating 
the Normal variates randomly, which we discuss briefly. Let 
be a Uniformly distributed random variate defined over the 
interval (0,l) and be a standard Normal variate. 

The first approach, discussed in Fishman [lO], uses the 
following formulae for generating X^. 

COS 271: 

X^ = (-2 logg sin gTcUg. 

The second approach, given in Naylor [ll] uses 


= ( ( 


k 

Z 

i=l 



X 
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n , " 

This is referred to as central limit approach . The number 
k, the number of to be generated, must be chosen depend- 

ing upon the required accuracy and computational efficiency. 

A high value of k is desirable, say around 24, is desirable 
for accurate results. But for computational simplicity, one 
may use k = 12. 

The third method, discussed in Laxv and Kelton [12] 
is as follows! 

Step Os Generate U-, and 

Say = 2U^ - 1 and ~ 1 

Set W = ^ y^2 

If (W > l) then go to Step 0, else go to Step 1, 

Stejp_ls Put Y = (-2 logg ¥)/W 

= Y and Y. 

The above three methods were programmed and histograms 
were plotted by generating 1000 random variates for a given 
parameters of Normal random variable. The shapes of the his- 
tograms obtained with the third method were observed to be 
very much similar to the Normal distributions of input para- 
meter. Hence, the third method was found to be most appropriate 
for generating Normal random variates and used in the computa- 
tional study. 
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(ii) No. of Simulation Runs in Simulat ion A-pproac hs 

For olDtaining the mean and standard deviation of 
project duration using simu].ation approach, the network was 
simulated for a total of 1000 runs. It was observed that the 
values of the parameters get stabilized after performing 
1000 Simula!: ed trials on the network. 

( i i i ) Generatio n pf_ th e^ . Par ame t e r s_ p I.. j - 'vI'^'V’ I.Qfls. s 

After generating the network structures using the 
network generators, the means and variances of durations of all 
the activities in the network are to be generated. These are 
generated in three ranges individually. The first range consi- 
dered is (5, 10), which is selected as narrow in order to 
clearly bring out the effects of near critical paths in the 
network. The next range is chosen as (l, 100 ) in order to 
make the situation closer to more practical situations. The 
third range considered is (5, 50), 


(iv) 


The V alues pf_ s and Limit 


is, P ath_ Enumeration Algorithm s 


In the Path enumeration algorithm, a small value of e 
is used to check for the stabilization of the parameters of 
project duration. The stabilization of these parameters depends 
upon the number of paths enumerated, which inturn depends upon 
many factors such as network, size and structure, parameters 
of individual activities etc. Thus the value of e is correlated 
to a given instance of the. problem. In an effort to relate e 
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with the problem, by a constant observation of several problems, 
it was found that the parameters of project duration would not 
change once the difference between successive values becomes 
nearly 0.2 percent of the mean of the first longest path for 
four or five times. Thus the value of e was fixed as 0.2 percent 
of the mean of the first longest path in the network and it 

was. checked, for termination, whether the difference between 

\ 

successive parameters is less than s for five consecutive times, 
i^so, the parameters were observed to stabilize after the ehu~ 
meration of only a few paths in the network and in most of the 
cases, before enumerating 40 paths. Thus the limit for the 
number of paths to be enumerated was put as 40. 

4.4 COMPUTATIONAL RESULTS! 

A computer program was v/ritten in FORTRAN- IV for the 
Path enumeration algorithm. The program stores the network 
in adjacency array representation. The conventional method, 
simulation method and Sculli 's method were also programmed. 

The listings of the programs are included in the Appendix. 

These programs were debugged and tested on DEC-1090 multi- 
programming, time-sharing computer system. 

The main emphasis in the computational study was Taid 
on (i) the accuracy of the results, and (ii) the computational 
times taken by various methods, Ifl/hile solving the Path enumera- 
tion algorithm, the number of paths enumerated for reaching 
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the results is also emphasized. Since these factors depend 
upon the network size and its structure, a large number of 
problems generated v/ith the three kinds of network generators 
NG 1, NG 2 and NG 3 were used in the study. The emphasized 
factors also depend upon the means and variances of individual 
activities in the network. Hence three ranges, as discussed 
earlier, were used for generating these parameters. If (kl, k.2) 
is the range used, the parameters were generated as follows* 

m.. = R/d© (kl, k2) V (i,j) e A 

X J 

V.. = RAI© (kl, k2) ¥ (i,j) e A 

X J 

where RAND (kl, k2) is a function that generates a uniformly 
distributed random number in the interval (kl, k2). 

The accuracy in the parameters obtained bj)^ each method 
is measured by taking the percentage deviations from the 
simulated values. Since for large networks, with a higher 
value of project duration, the percentage, deviations are not 
appropriate for measuring the accuracy, the actual values 
obtained by each method are also shown for each problem. 

To generate the network problems by NG 1, the length 
of the network v/as varied from 5 to 30 and v^idth from 5 to 8, 
All these problems were first solved by fixing the range for 
the parameters of activity durations as (5, 10). Table 4.1 
shows the results of these problems for various methods. The 
same problems v/ere then solved for other ranges (5, 50 ) and 
(l,100). These results are presented in Tables 4, 2 and 4,3 
respectively. 



48 


In the network generator NG 2, the parameter n was 
varied from 5 to 60 and the problems generated were solved for 
all the three ranges (5, 10'), (5, 50) and (l, 100), These 
resiilts are presented in Tables 4.4 to 4,6 respectively, 

A number of problems were generated by NG 3 also by 
varying the width from 3 to 10 and length from 5 to 30. The 
number of arcs were chosen reasonably to relate the problems 
with real life problems. Tables 4.7 to 4.9 present the results 
for these problems. 

To have a clear picture of the errors in the parameters 
with various methods, the percentage deviations are plotted 
for each problem. Fig. 4.5 shows these deviations for the 
first six problems shoiAmi in Table 4,1, The errors in mean 
are shown in Fig, 4.5(4) and the errors in standard deviation 
are shown in Fig. 4.5(B). Similar kind of plots are shown 
in Fig, 4,6 and Fig, 4.7 for the network problems generated 
by NG 2 and NG 3 respectively. 

The execution times for various methods are noted as 
f ollov/s, (the execution time is the CPU time exclusive of input 
and output times )s A network is generated with given dimensions. 
Five different problems, for the same network structure, v/ere 
generated by varying the seed from 0 to 40 in steps of 10. 

These five problems were then solved and their average execu- 
tion time was noted. Tables 4,10 to 4.12 shows these times for 
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the network pi'oblems generated by MG 1, NG 2 and NG 3 respec- 
tively. The average execution times of both conventional 
method and Sculli 's method are alv/ays lesser than those of 
Path enumeration algorithm. The results shown in Table 4,12 
(for NG 3) are plotted and shown in Fig. 4,8. 

The computational experiments prove that the conven- 
tional method, as anticipated, always under- estimates the 
mean and over-estimate the standard deviation of project dura- 
tion, The errors are more when the parameters of activity 
durations are in the range (5, 10), This is due to the effect of 
many near critical paths in the problems generated with that 
range. Sculli 's method always over-estimate the mean and under- 
estimate the standard deviation. This was also expected because 
the dependency among paths was not considered in this method. 

The errors are more with the dense networks generated by NG 2 
because of more number of dependent paths terminating at each 
node* 

Path enumeration algorithm performs as a compromise ' 
between the above two methods. The deviations in the parameters 
may fall on either side of simulated values. The reason is as 
follows s The algorithm starts with the results of conventional 
method and improves them gradually by taking the effect of near 
critical paths, until the successive values of the parameters 
differ by a negligible amount. In this process, the mean of 
the project duration ma^^ exceed the simulated value by a small 
amount. 
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The improved results with the Path enumeration algorithm 
are obtained at a cost of extra computing effort. However, the 
total effort is still v/ell belov/ the effort required in simula- 
tion method. The algorithm requires a higher computer storage 
when compared with the storage requirements of other methods. 

This is due to the reason that Dreyfus method requires to 
store all the previous paths when it evaluates k-th longest 
path. This limits us to solve still larger size problems. The 
computer program occupies mk + 3nk + 8m + 3n words of central 
memory. 

4.5 CONCLUSIONS; 

In this thesis work, a new approximate method for deter- 
mining the completion times of PERT networks is developed. Com- 
putational analysis over a large nimiber of randomly generated 
network problems shows that the method will prove of assistance 
when the accuracy in the parameters of project duration is 
important. The method can solve practically large problems in 
reasonable amount of computer time. Thus it may be concluded 
that the method is 'appreciably-acceptable' for estimating the 
project duration in PERT networks. 

Though it is assumed that the activity durations are 
Normally distributed, the method can be used for any kind of 
distributions, even with different distributions for different 
activities, because the paths will always be Normally distributed. 
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TABLE 4.10; Execution times for PEA and Simulation (for KG l). 
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'^RETORii 
■ EMO, 


•: rr:} 

REA3U3,*), ■1ii:c-l:l'/,.Si)COviV,ij:;SI^;a,v^;JSI '..viUr^J, . SOHjRt , 

1 ^c;3cai,,3i).sci.'u 

RlC3M/= ir:2j',i V/ ' 

R2C0'i\/=S>J20;i^/3i)3l;lU 
R1 HurUs^iEHUR l/il3Si ii) 

R2HJR.l=5l)HlTR i/si>3I ''iU 
R1SC3Ci=>1ESC0L/mh..SI-i!} 

R2SCJl4=SDSCUL'/3i)3Xn!J 

WRI rs ( 22,10 j , 1 CUNV , H/ICO'^ V , lUHOra , K2:-l uiu » 3 1 3CU L , H2 iZ 31. 

R I P C = h.ni ( 1 . ' > - 1 C J 3 i' ) X I 0 . ■;, 

R2PC=ABS ( I v<> -p'2 CO': V ) *. 1 j'? . C ■ 

Hlpt-i?AB3( ,l.O-Rl HURD^IO'J.O 
R2P»i=ABSCl.B-R2rUIRl)-*lOj. J 
R:1PS = :VBS{1 ,0-Rl 3CUli)>10^,u . , 

R2p3 = A3S(l.t»-ft2 3C‘J[4)*lO':/. ') 

vIRl: i’S C 22 > 15), R1 PC> r!^PC , Ri PH , R2PH , 2,1 PS >.H?PS 

roR'Mra5^X»6fB:.2.)^' 

.aEfURi.' ■■ 



INTElSeR Aans ,TAlLt90.a ) , hTAD ( 9yv ) , PDI A‘T( 30 0 J r( 301' ) , > 

1 p.Lisrc9r>o) / 

REAL ME’A:'<'f.':K'u),VAK(i<0'>3 

c: 0 M n A. 0 • VE S , A R c 5 ,> r ; 1 0 f ' i fc A D , P u I , , i-' P D I U T , K I., J. S T , £ A ■ i , V A i? , 1 y E 
CALI. :'Gl . 

PAUL S3 M J 

CA&L PEA • . 

CALL SCULU , ' ' 

DPE^U ’IP iTsi;3 , 'MJPKs ' ASC 1-T: M 
CALL PATIO 

CLO L E ( t)M I*T=2 J , M noCs/ AvSC U * ) ' 

STOP 
BRO ' ' 

THIS SPoK-OilTtOE GE'LERATrS' AKTKnpKS OF TTPE;-! . 

.SUS'ROUrTRE ,OGT- 

I « r E GE R ’ i’ f-.n T K , AR-C S , C 0 U N T, , T A1 L ( 9 0 0 ) , H E AO ( 9 0 0 ) P Oi R V C 3 0 P ). ^ 

1 AL.lSrC30M!) , Bl4ST( 300 ) rHPQINtC 30o ) o 

■ real. 4EAH-(90>r),VAR( 4)0.;);'- ■ . a'' ^ ''’’a; O''' : 

CDH9 3i'i r/aOES » ARCS , TAJ L , HE AD rPOTPX, PPO I AT f RLI .^T , MB A>K>. VAP , IS:B)5;l|: 
TY.P-E. I'O , ,'V 

rORyiATdOX# 'EOTFR THE tvIDlH REC>(.'T«Ei) •: ' ^$5 * 

ACCEPT "*VOin'TH . 

. TYPE'. 20' ^ ' .0'. " , , , 

FORBATt lOK, 'E hTER THE LENGTH REOHIREO : ' » «) 

' .ACCEPT <‘’>LE0.G'TH '■ . 

CALL SETHA-iCISESi;)) 

M4DES=LE;MGXH*wIOTriT7. 

" •couMr=L.-'.o 

Ob i'O As2,>aprH + l^ 

COU-OfsCOOilT + l 
TAlL(CPU«rJ=l 
HEA,bCCbua'.t) = I." ; 

.,;CONTiailE'.:, ,::o 
' ; -p0':;.''5O ^ ■KsTy'L'Ei’JG.T'H*'!..-;- 

;;tPa.rE;:g4T>G^,^0v^ 

■Vc0bB;rRppwr+i; 



HEAOtCLJif ir) = t4i 
COUf^rsCDHtiT+.J 
l’Alt.(CDd iT) = I 
HF,RvV(cni.)'iT)si+.. rorn 
IFd.Fv'.lDGu T) 5v 

coo.vrsc.iM Tfi 
TAIL.(C'M i l ) = I 
HE A 0 ( C 30 1 f ) = I + I u f rl f 1 
C.0MI:IU£ . 

IS»!lS + >iXr' Td 
irsir+wiDTd : 
pc 7 Cj I = ISri r 
IFd .EP,! OOP T3 6t/ 

CC>UNraCPU:iT4l 

xAitiCcoupnai , : 

MEAD(CDUdn = I + l 
CCPM’sCP'd^i'r+l ; ; ■ 

TAl.L(C3U;n'i = I 
HEAD(CPd:M nsi'lPDKS 
‘COHTn'tlf? 
ftRCS=CPu.';r 
DO 1=1, POPES 
RLisrcipse * 
coNri'jiJE 
poiPid j*i. 

'K=l' 

do: 90 J=1,AHCS 

ALlSr cheap (P ) ) = ALTSn iiEAD r JO ) 1 1 
IFCl‘ArU(J).ep..K)GO TO 90 
.■X=K:fl ■ 

' .'pOl.Pi’CK)*!, ; " 
r-coPtiNoe'' ' 

PClP.r(wPOES) = ARCS + l 
pCtIPr(MOOES + l) = ARCS-H 

apisrci) = 1 

\RPpiMi’d')..'ad 

^:po: P-0.C-I-®'2'^O0pfcS;'-; 

B LI s rd psiSP I s-rt i - d + al i stc ^ 

'■■RpOt^MXCdsBLISTll;).:' 

f Rpo t^M H i>OS$d ) « ARC S + 1 
pa 1 14 V j;alf ABCS 
■>**£ A^: t (5 <5 -liBi; 



via 

uo 

C- ■ 


BLlSrCMSAirCO J) = 3LIS'1'( J>) + 1 ; 

COwili'HJc; , , . 

irm TJ: E ( 2 2 , 1 2 ;; -) r ' i 1 ■'» T H , I 4 e/. • c 7 , .f-< 0 U F S hV ' < C 3 

FORHATCrOX, *W10 J’il = ' , T i , I E?:GTti = ' , I 3 , ' t):)3SS = ',,T3,- 

1 ' ARCS =- 

RFTUR^J 

EMD ■ ^ 


FUl’J-VlO '! FnP. GF'N!ERhTi;:J(J A lytS rRlBnxED R A >''aM ;’ JABER 

Bei‘"'’EEv TviLy'SPBCIFiSD Xi'iTBGLRS. : 

FUfjCTlO.-' RAfjDCtCflC') 

IRA^^=I: + (^r;-tC•^i ) *RAri(X) 

IFCIRA' ,G ",ia)IRAvi = 10 . 

RAry.[) = IrvA!J ■ ■ 

RFi-^Rfr . 

■ Ef-jy ■ 



A«CSi!TAIu(t60'’),Kg;AD(16Ov),V3X';T(60J ,KP3I.'Jr(50), 

I: RL.ISTi[15,0G)' , 

RRAU .''lE'AN(ir>o'.-).,VAHf IfSOf ) 

'Ifc; A 'J , VAR , ISEED 

C/itli*' '•I'J? 

CALL’ SlAO 
CALL' PEA 
CALL, 'scuun-. 

OPEVtU-xil fs21, .ia.f)E=:'AS.CII ') 

■CAUV-RArn ' . ■ 

qLOSE( J^aTa23,MaijesVASCII'J^ ^ 

.^SIDP’ ■■ 

,.EV0 


THIS S^BROUTlRe rGEf^ERATRR R51V.'0»K;i, JF TYPE-2. : 

.suBROirriNS.' --iGa 

INTtSER Yi;tas, TAli,(l t,uca vBEAD(A6<v>) ,c y>i.Vr, ALtST(> j ) , ALtSTCoO ) > 
1 pap i'( olV) , RPOI W Ti 60 ) » RUST (1600) 

. REAL ■ n S A^■l (,T 6 y AR ( 1 60 n.) ■ ^ 

^oaO£SrABCS., ^AIli|^^e:AD# PDl :vT,RP;Jly iSRuIST/'IEL^ 

, TYPE' 10 ■'■ ' 

. FOR.HAT(iaXpe:a£R THE of rvOOES : ' ,s) 

accept ..♦,TU bEvS ■' , 

■co'.i:'jr=o', - . V ■ 

■.JJsSlO'DCS' ■' 

DO 20 Tsl,,-i-t 
00 20 

CQii:'ir*CpUMT+l 
rAlLTcaURDsI 
(IKAOCCOlIfaXsK 
Cor4TiriUE . 

■ARCS:«'C0;.1NT'..':.''; 

■ 0:0'' -30' 'T^iivLaaE ." 

^Afrsrcrpo' 

CO 'r' I V'V''2 

,'pO:X'N'T(l):=i:; 



;'.4'0 COM.riNOE\,' 

pdi<'ir(:^oDe.s)saHCSfr 
poi >r(iNiU(>ssf i)3kf'Cs+i 

0 L,isrti) = i 

f^poiNixd )=i 
DO 5u. I=?,:<'juc;3 

BLiSr ( I D^sbUi. ST t J:-l )+h(ilS r(J-l > 

: 50 ■ , CONrjriOE" 

RPOIi'>JTCjv'OOE:Sfl )sAKCK + l 
00 60 J = 1,AP.CS 

tfAR( J) = k4 !D(5,10) 
«iiJ.5C(4t,rsrcbb:A!>(j) ))=j 
6LI s r ( H EAO ( J )■ D =4Ia 5T ( H?:ALHd ) ) 4 1 
60 ■ ' ca'iii.^niE ' ' ^ ■ - ; 

£ r S,( 2 ? , 7 0 ) /ri ij D ss',: f. RC S 
70 >OK'14r(lOiC,'>UiOt5 = ',13,' ARCS = 
RErUR?i' 



liiTSSeH A'^Cfi ,ThIL( 1200) , HEAOC 1 200 j , I’DI’J T ( 1 00 ) , SP'^I '*TC 1 00 ) , ' 

■ 1 1 o(..isr,ci2oo J ' 

■ - ^AL 

COM'O j r/aL,r,EHn,POXviT,OOH.H’,f<LIST, i£A J,^A0, ISESD 

.' z\hi o:;? , ' ■ ■ , 

. CALiO' ai:'1U . ' ' ' . , ■ 

;C.\bf4. or:.v ■, ■ 

CAfjO 3C0t.f.I , ' 

OPE ') C I -i t r=2 i /■'■.OOEs ' ASC It') 

CALL. .lATIO - 

•CLOSeCO HTa^ J, vl'JOKa'AiiClI') 

'sro*^, 

■ SNIO ;.. 

c T‘1TS S-.MSOiJTI.'JE v^eUEftATES rrETviOHK.s OF TCPS-J, 

; SaOROUTMS,' OG'o", ' ■ ' ■ , ^ ' 

IlTECEL ARCSf AlOrfl,0'AIL(12(»O),HSAc0(l2OO) rPOintC lObOffiPOtlTt 
: \ «i4 13 r (1200 ) OO ’ir , ACiOO ^ lOO) r AL13T ( 100 ) ,}0H3r( I 00) 

■ ';iSAu ^«A.i(1200),0AR:(l-2x)o-) v',' V ' 'v,;- -V-'; 

COIlMO’J vdo^:S, ARCS/ TAIL, HEAfr, POlOT/KPDlfi r, RLt Sir '^SA'irVAR : 

■ • 'type 'ro' 

10 PORRA r (1 OXr 'EriXER THE WIDTH ! ' > $) 

ACCEPT '' . : 

: type- '2C' ’ ■ ' : ' ' ■. 

20 POROATlIL^, 'E'- rfEE the LEHETrl. ;. V,$) - 
accept * /LS*!3TU, ’ 

.. 'ilfl = ((,EN^3TH + 10^/Hi)rH 

^ f^ODES=v.'IOTH»uEE3TH+2 
, .■iAK=(O0D'E3* C.'*'0i-'E'S-) 0;)/2' ' 

.type. '30r:4lH,oAX ' 

30 FnRVlAi’(l0X:r*E0TSR TflE iiO.OF ARCvS t?i titE RAiMGE" , 13, ‘ - ' f ti) 

\ .TYPE' 'it'V '.-V'. 

'40'. , ■■P'0RR4T(7>'1 C.'Xr*Hi'iTEK:;THE-,i'U/::OF ■-^ARG3^':■^ %3) : 

'"■'"accept..-'*:, AR'CS".' ;■/> Ax".. 

..'COURT-^'.'v^" 

■'^A.’:n7..':.obA5o^;:j=2,ViOT 
:' . 'x' 'x''bbU'i^raCOUWT.+rlX''-;;;^.;-' 

v''sb';A'cbbT.i.?i:oEAv 



, , 00 60 , JsK^K'* .7 CDT.-I-I 

cnu.')r=C'i(j^r+i V; ^ ; 

A( Jf J + = l 

_.60 . CDfirXi'^ir. ■ , ' ■ : 

: : KsK+i^ior;! 

■'TO',-,' 

DO 'i>0 . rs,^yK + .ili>i’/>l 
ca(JN.p=CM:jiv+i' . 

Ad 

■ ' 80 '■■ . co^idiiutn ^ 

. 90 IFCOOnT.SR. ARC3J GOTO 140 . 

100 'NDOElrHAiilOCl^rJOOGS) 
llO NOOF2arvA;iO(l /NO.UKS) . 

IF't 'dOcOJ - iv)pe2) l?Ud'l^d36 

120 iF(Ati'ir)Dt:i » ;i 01)152 ).kcu1)gutov'^'^ 

' ■/,. ■ cjU''-)r=oj(i':JiT+i ■' /■ ; 

■ AO'iOOFl ,;i3aF:2) = l V - 

G0T3 'Of/ ■ ■ 

J 30 ^ IF( A CMOOe^ , JODKi ) . sOQ. 1 > GOTO 100 ; 

' 1, 'coo>'J'r=vair;iT,+ i. 

A-(N0l)S'2,:i:jDFl;) = r - . ' 

GOTO; '9 0 ; ■ ■ 

' 140 AOa':J raO 

DO 150 1=1 , MOORS 
OD ISOJSJ, MOOES : 

IFd.GG.J) GOTO I5y 
IP(A(r/j).ME.l) GOTO 15(; 

IC00'5i' = K00MT+l 

' > 

'd'Alu'dCOUM.Tlst ' 
ilSAOCXOUMds'J 
■:150'',: COMj’lMO'E- f v. 

If (XOU aI'.EO.APCS) GflTO i'7U; 
type IbC- , ■ ' 

161E f^Lip^ATCIOX, 'SOME /JRONG IN TuE PROGRAM. 
' ■ STOP';.; 

170 DO Iso I: = 1»?'0DES ; 

, 'A'LlS'dl')at'': , 

■ '■ i8'0 . ;■ ,V 

V.\V' 

^5. ;pGXNr'tl)al' 

' ''do.',' 1'90' 7-i ,''ar.c;s,'.-'' ■ '.■ 

V'' ■'^'^;^LlSY'('rtE,AOt;0;3 l=AW'Sf (li'EAD.ti'j^ 
■';iFC'rAxG(a^Eg;K7:'GO’;?TO;:i90,;,;;;v^ 

MlB 


190 co^:M. , 

; POlNtC I.lOKSJcARCS+l 
pai’'.:r( it)ne:s+i|s/^vRC5,f 1 
Bhisr'C 1 )=,,!■ ■■ 

' 'RP3I'ViT(i )s) 

• 03 T = ?#'.:j'>i:s 

Bi.rsrt f )svi[,i3ri T-i 3 + AL,tsTu-i) 

RpOI.^rCDsLiUlSTU) - 

200 . cnwi’iriu'f: ' , ' ^ 

Rpp) *iT(,!OOt:S + n=AKC.S+l 

■D0.:2if- jsUARCi;. .. ■■ ■..■:'■■■■ 

■ MEA.N(J3 = HAwp(5,10}^ _ 

yAfaj)sKA'io(5,,i to 

■W ■ \ R’wISr(8{,ISTCHErtO(JJ))aJ 

Bt/ISJ'C4t':/U)( J ))=:biaST( BK/U'I J j ) + l 
■. 210. ■ CCibl'I'-ibiv ■ 

WRire(22, 2'2oy^wl!)TH,tif:Kn'rH,':500F:S, ARCS 
220 FaBN.AT(iCX, LEfiGXn » HOOfiS 

1‘ I4»' ARCS = ',15///) 

PEtliRf] 



